1 Pendulum 



1.1 Free oscillator 

To introduce dynamical systems, we begin with one of the simplest: a free 
oscillator. Specifically, we consider an unforced, undamped pendulum. 




The arc length (displacement) between the pendulum's current position and 
rest position (0 = 0) is 

s = W 

Therefore 

s = W 
s = W 

From Newton's 2nd law, 

F = mie 

The restoring force is given by —mgsind. (It acts in the direction opposite 
to sgn(^)). Thus 

F = mW = —mg sin 9 

or 

— r + -sme' = 0. 
dt'^ I 

Our pendulum equation is idealized: it assumes, e.g., a point mass, a rigid 
geometry, and most importantly, no friction. 
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The equation is nonlinear, because of the sin^ term. Thus the equation is 
not easily solved. 

However for small ^ <C 1 we have sin^ ~ 9. Then 



whose solution is 



g 



or 



9 = 9o cos{iut + 0) 
where the angular frequency is 



the period is 



T = 27rW-, 
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and ^0 and come from the initial conditions. 
Note that the motion is exactly periodic. 

Furthermore, the period T is independent of the amplitude ^o- 




1.2 Global view of dynamics 



What do we need to know to completely describe the instantaneous state of 
the pendulum? 

d^ 

The position 9 and the velocity — = 9. 

dt 



Instead of integrating our o.d.e. for the pendulum, we seek a representation 
of the solution in the plane of 9 and 9. 

Because the solution is periodic, we know that the resulting trajectory must 
be closed: 

e 




e 



In which direction is the flow? 
What shape does the curve take? 

To calculate the curve, we note that it should be characterized by constant 
energy, since no energy is input to the system (it is not driven) and none is 
dissipated (there is no friction). 

Therefore we compute the energy E{0^ 9), and expect the trajectories to be 
curves of E(9,9) = const. 



1.3 Energy in the plane pendulum 




The pendulum's height above its rest position is h = I — I cos 9. 
As before, s = arc length = W. 
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The kinetic energy T is 



2 2 ^ ^ 2 

The potential energy U is 

U = mgh = mg{l — I cos 6) 

= mgl^l - cos 6) 

Therefore the energy E{0,0) is 

E{e, 9) = + mgl{l - cos 9) 

We check that E{9,9) is a constant of motion by calculating its time deriva- 
tive: 

dE 1 



= -ml\299) + mgl9sm9 



= ml^9(9^jsm9 
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= (since the pend. eqn. ^ = — -sin^) 

L 



So what do these curves look like? 

Take to be the highest point of motion. 




Then 
and 



^(^o) = 



E{9o,9 I^J = mgl{l - cos6'o) 



Since cos^ = 1 - 2 sin2(^/2), 

^(^0, \eo) = '^mgl sin 



2 / ^0 



= E{9, 9) in general, since E is conserved 
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Now write T = E -U: 



^mfe^ = 2mgl (^sin^ ^ - sin^ ^ 



For small such that 6 <^ 1, 



I 



l\4 4 



(1) 
(2) 



or 



e 



+ 0' 



e 



Thus for small 9 the curves are circles of radius 6*0 in the plane of 9 and 

e 




What about large? 
Consider the case 9o = ir. 

For 9q = IT, E = 2'mgl, and equation (2) gives 



9' = 4^ 

L 



Thus for 9o = IT, the curves are the cosines 



'9' 
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^-1012 

e/7c 

Intuitively, we recognize that this curve separates 

oscillatory motion [E < 2mgl) from rotary motion [E > 2mgl). 



Thus for undampled, nonlinear pendulum we can construct the following 
phase portrait: 




The portrait is periodic. 

The points 6 = 0, 6 = . . . , — 27r, 0, 27r, . . . are stable equilibrium, or fixed, 
points (actually, marginally stable). 

The points 6 = 0,6 = ..., — Stt, — tt, tt, Stt . . . are unstable fixed points. 
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The trajectories appear to cross, but they do not. Why not? 
(Deterministic trajectories.) 

If the trajectories actuahy arrive to these crossing points, then what happens? 
(The motion stops, awaiting instabihty. But we shall see that it would take 
infinite time to arrive at these points.) 
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